Stability of atomic clocks based on entangled atoms 
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We analyze the effect of realistic noise sources for an atomic clock consisting of a local oscillator 
that is actively locked to a spin-squeezed (entangled) ensemble of -/V atoms. We show that the use 
of entangled states can lead to an improvement of the long-term stability of the clock when the 
measurement is limited by decoherence associated with instability of the local oscillator combined 
with fluctuations in the atomic ensemble's Bloch vector. Atomic states with a moderate degree of 
entanglement yield the maximal clock stability, resulting in an improvement that scales as N 1 ^ 6 
compared to the atomic shot noise level. 



Quantum entanglement is the basis for many of the 
roposed applications of quantum information science 
l|. The experimental implementation of these ideas is 
challenging since entangled states are easily destroyed by 
decoherence. To evaluate the potential usefulness of en- 
tanglement it is therefore essential to include a realistic 
description of noise in experiments of interest. Although 
decoherence is commonly analyzed in the context of sim- 
ple models Q , practical sources of noise often possess a 
non-trivial frequency spectrum, and enter through a va- 
riety of different physical processes. In this Letter, we 
analyze the effect of realistic decoherence processes and 
noise sources in an atomic clock that is actively locked 
to a spin-squeezed (entangled) ensemble of atoms. 

The performance of an atomic clock can be charac- 
terized by its frequency accuracy and stability. Accu- 
racy refers to the frequency offset from the ideal value, 
whereas stability describes the fluctuations around, and 
drift away from the average frequency. To improve the 
long-term clock stability, it has been suggested to use 
entangled atomic ensembles [1 El an d in this let- 
ter we analyze such proposals in the presence of real- 
istic decoherence and noise. In practice, an atomic clock 
operates by locking the frequency of a local oscillator 
(L.O.) to the transition frequency between two levels in 
an atom. This locking is achieved by a spectroscopic mea- 
surement determining the L.O. frequency offset 6u> from 
the atomic resonance, followed by a feedback mechanism 
which steers the L.O. frequency so as to null the mean 
frequency offset. The problem of frequency control thus 
combines elements of quantum parameter estimation the- 
ory and control of stochastic systems via feedback 0, Q • 
The spectroscopic measurement of the atomic tran- 
sition frequency is typically achieved through Ramsey 
spectroscopy gj, in which the atoms are illuminated by 
two short, near-resonant pulses from the local oscillator, 
separated by a long period of free evolution, referred to 
as the Ramsey time T. During the free evolution the 
atomic state and the L.O. acquire a relative phase differ- 
ence 5(f) = SujT, which is subsequently determined by a 
projection measurement. If a long time T is used, then 
Ramsey spectroscopy provides a very sensitive measure- 



ment of the L.O. frequency offset 5ui 9]. Here, we inves- 
tigate the situation relevant to trapped particles, such 
as atoms in an optical lattice or trapped ions 
In this situation, the optimal value of T is determined 
by atomic decoherence (caused by imperfections in the 
experimental setup) which therefore determines the ulti- 
mate performance of the clock. 

We consider an ensemble of N two-level particles with 
lower (upper) state | J.) (| |)). Adopting the nomencla- 
ture of spin-1/2 particles, we introduce the total angu- 
lar momentum (i.e., Bloch vector) J = X^=i Sji where 
e.g. Si = (| T)j<T I - I I) A Initially the state 

of the atoms has mean (J) along the z direction and 
(Jx) — (Jy) = 0. Unavoidable fluctuations in the x and 
y components (J 2 ) = (J 2 ) ^ 0, result in the so-called 
atomic projection noise. These fluctuations give rise to 
an uncertainty in the Ramsey phase 5<t>p ~ AJ y /\(J z )\ 
as indicated geometrically in Fig. 1 0, [u| . For uncor- 
rected spins aligned along the z axis, the uncertainty 
from independent spins are added in quadrature, result- 
ing in the projection noise AJ y = y/N/2. To reduce 
the measurement error it has been proposed 0, ll'l FLU 
and demonstrated 0] to use entangled atomic states (so 
called spin squeezed states), which have reduced noise 
in one of the transverse spin components (e.g., J y ) and 
non-zero noise AJ Z in the mean spin direction. Ideally 
this gives an improvement by a factor £ = ^f~N AJ y / (J z ), 
which can be as low as £ = 1/VA for maximally entan- 
gled states 

Using a simple noise model it was shown in Ref. 
that entanglement provides little gain in spectroscopic 
sensitivity in the presence of atomic decoherence. In 
essence, random fluctuations in the phase of the atomic 
coherence cause a rapid smearing of the error contour 
in Fig. la. For example, dephasing of individual parti- 
cles results in an additional contribution (N / A)(8(f) 2 ) to 
the noise, where (5(f) 2 ) denotes the variance of the phase 
fluctuations (increasing with T as (5<fi 2 ) = 7T for white 
noise, where 7 is the dephasing rate). In practice, the 
stability of atomic clocks is often limited primarily by 
fluctuations of the L.O. As we show below, the L.O. flue- 
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tuations result in the added noise A J z (5(f) 2 ), where A J 2 
is the initial variance in J z . This added noise is due to 
the error in the feedback loop, caused by the longitudi- 
nal noise AJ Z . For weakly entangled states, the added 
noise is considerably smaller than in the case of atomic 
dephasing and the use of entangled states can lead to a 
significant improvement in clock stability. 

In what follows we outline a model that incorporates 
the effects of atomic noise and spin squeezing as well 
as that of the feedback loop. Before proceeding, we 
note that qualitative considerations along these lines were 
noted in Ref. 01- At the operating point, the error sig- 
nal in Ramsey spectroscopy || measured at time t k is 
determined by the operator 

N 

E{t k ) ~ sin [<%(*fc)l + % cos[c%(t fe )], (1) 

3=1 

where 5(f) j (if.) is the phase acquired by the jth atom dur- 
ing the interrogation time T and all operators refer to 
the initial atomic state. We separate the phase into two 
parts 5<fij = 5<f>o + 5<jP E , where 5(j>o(t k ) = Jq Su(t k — t)dt 
is the phase due to the frequency fluctuations Su> (t) of 
the L.O., and 6<j) E is a phase induced by the interaction 
of the jth atom with the environment. In order to lock 
the L.O. to the atomic frequency, the interrogation time 
should be short enough that (S(f>(t k ) 2 ) < 1< Expanding 
in terms of S(f>(t k ), we find the measured error signal 



£(t h ) ~ (J z ) [6<t>o(tk) 
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N 



/ 3=1 

+ [Jy(tk) + SJ z (t k )Wo{tk)} + ■ ■ ■ ■ (2) 

Here SJ z (t k ) = J z (t k ) - (J z ), where J z {t k ) and J v (t k ) 
are random numbers with a distribution corresponding to 
the initial atomic state (we consider here states for which 
{JyJz) = 0, so that we may treat J z {t k ) and J y {t k ) as 
independent random variables). The term multiplying 
(J z ) in J5J) is used to estimate the frequency offset, while 
the remaining terms represent measurement noise. 

The feedback is started at t = and, at the end of 
each Ramsey cycle, at t k — kT (k — 1,2, ...), the detec- 
tion signal is used to steer the frequency of the oscillator 
to correct for the fluctuations accumulated during the 
last cycle 6uj(t~£) = Suj(t^) + Au(t k ), where and f£ 
refer to before and after the correction, and Auj(t k ) is 
the frequency correction. Assuming that negligible time 
is spent performing the 7r/2 pulses and in preparing and 
detecting the state of the atoms, the mean frequency off- 
set after running for a period t = nT is then 
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We begin by analyzing the simplest case of linear feed- 
back (in £(tk)) and later extend to the more optimal 
nonlinear feedback case. With Auj(t k ) = , using 

(|2~|) and substituting in 0, we find, ignoring for now the 

6(f>o(t k ) 3 term, 



5Q{t) 



-1 ™ 

—j- YWvQk) + 5J z (t k )Sq>o(tk) 
T \ J *) k=l 



N 



3 = 1 



(4) 



Note that the acquired offsets S(f>o(t k )/T (k = l,...,n) 
due to L.O. frequency fluctuations are corrected by the 
feedback loop and do not appear in while measure- 
ment noise is added at the detection times t k . The first 
two terms in (@J are uncorrelated for different t k since 
the atomic noise for different detection events is uncorre- 
lated. If the dephasing noise is uncorrelated for different 
t k , then the fractional frequency fluctuation (Allan devi- 
ation) a y (r) = ((^(r)/^) 2 ) 1 / 2 , is 



AJ 2 2 W 2 ) 



(\(J 2 z ))(Hl)] 
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Here A accounts for the possibility of collective decoher- 
ence, so that for atoms dephasing collectively (indepen- 
dently) A — > 1 (A — > [N/A)/(J 2 Z )). The L.O. noise affects 
the atoms in a similar fashion than collective dephasing. 
Note, however, the significant difference between collec- 
tive environmental dephasing, which enters expression 
© as (J Z )(S4> 2 E ), and L.O. noise, for which AJ 2 ((5</>q) 
is the relevant expression. The feedback loop results in 
a large cancellation of the effect of the L.O. noise on the 
stability; the uncanceled part of the noise is now propor- 
tional to A J 2 < (J 2 ). 

When decoherence is negligible, (S^q) = (5<j) E ) = 0, 
the long term frequency stability is given by cr v (r) — 
AJy/ujy/rTiJz) as shown in Refs. 0,0- F° r an uncorre- 
lated atomic state, the stability improves with increasing 
number of atoms as TV -1 / 2 10] . The maximum possi- 
ble improvement using spin-squeezed states is a factor of 
N- 1 ' 2 , yielding a stability <t v (t) oc N _1 0]. 

The best long term stability is obtained with the 
longest possible interrogation time T. When the inter- 
rogation time is limited by environmental decoherence, 
the latter cannot be ignored. This corresponds to the 
situation considered in Refs. [lfil EfT ]. in which case no 
substantial improvement is possible. In the practically 
relevant case where the main source of noise is from the 

l.o. muni 

the situation is quite different. In this 
case it is undesirable to use a very highly squeezed state 
with A J y ~ 1 because it has a very large uncertainty in 
the z-component of the spin AJ Z ~ N, which accord- 
ing to Eq. © has a large contribution to the noise. A 
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FIG. 1: a) Representation of the probability distribution on 
the Bloch sphere for a spin squeezed state \tp(K)), with k, — 
7V 1//4 corresponding to the squeezing parameter £ = jV _1//4 
(TV — 10, both the initial state and the state just before de- 
tection are shown for clarity). Thick lines indicate initial un- 
certainties in J. b) Noise spectra due to L.O. frequency fluc- 
tuations when free running, when stabilized to unsqueezed 
atoms (£ = 1), and when stabilized to spin squeezed atoms 
(£ = iV~ 1/6 ), TV = 10 3 and 7T = 10~ 2 . c) Inverse frac- 
tional frequency stability l/a y (arbitrary units) vs. Ram- 
sey time for white L.O. noise, TV = 10 J , with linear feed- 
back to uncorrelated atoms (o); linear feedback to correlated 
atoms (o, £ = 7V~ 1//4 ); nonlinear feedback to uncorrelated 
atoms (•); and nonlinear feedback to correlated atoms (filled 
o, £ = 7V -1//4 ). Points: numerical simulations, lines: analyti- 
cal results. 



moderately squeezed state can, however, lead to a con- 
siderable improvement in the stability. This observation 
is the main result of the present Letter. 

To find the optimal stability, we first optimize (jSJ) with 
respect to the interrogation time. Considering uncorre- 
lated atoms first, we have AJ y = \/~N /2 and AJ Z = 0; 
Eq. J5) then predicts that a y (r) decreases indefinitely as 
1/y/T. To derive Eq. however, we have linearized 
the expression in Eq. ft). an d this linearization breaks 
down when the (neglected) cubic term in is com- 
parable to the noise term that we retained, i.e., when 
5(/)(tk) 3 ~ AJ y /(J z ). In a more careful analysis 
based on Eq. including perturbatively the nonlin- 
ear terms in a stochastic differential equation, we find 
the optimal time 7T = (2AJ 2 /( J z ) 2 ) 1/3 . At this point 



C= — A 1 / 3 

<> 2 4 /3 




1/2 



(6) 



To evaluate the potential improvement in stability 
by using squeezed states (i.e., the scaling with increas- 
ing number of atoms A, in the limit A 3> 1), it is 
convenient to use a family of states parametrized by 
a small number of parameters. A one-parameter fam- 
ily of states that includes the uncorrelated state as well 
as spin squeezed states is given by the Gaussian states 
\iP(k)) = A^K)£ m (-i)' me ~ ( ' m/K)2 |™), where \m) are 
eigenstates of the J y operator with eigenvalue m and the 
total angular momentum quantum number is J = A/2, 
and jV(k) is a normalization factor. The transverse noise 
for these states is given by AJ y = k/2. For a large 
number of atoms A 3> 1, the uncorrelated state is well 
approximated by \iP(k = y/N)), while highly-squeezed 
states are obtained when n — ► 1. Within this fam- 
ily of states the optimal value is C — L42A r_1 / 6 for 
k ~ 2 1 / 16 A 1 / 4 (£ ^ N- 1 ^) giving a stability scaling 
as A -1 / 2 . This represents an improvement by a fac- 
tor of A 1 / 6 compared to uncorrelated states, for which 
C = 3/2 4 / 3 and the stability scales as A -1 / 3 . We em- 
phasize that these results are derived assuming a linear 
feedback loop. 

To confirm these predictions, we have made exten- 
sive numerical simulations of the frequency control loop, 
along the lines of Ref. @|. The noise spectrum of the 
free-running oscillator is defined by S(f)S(f + /') = 
(Su)(f)6u)(f)), where Su)(f) is the Fourier transform of 
the stochastic process Sui(t). We generate the corre- 
sponding time-series and at the detection times t/. = fcT, 
the accumulated phase 8<f)o{tk) is calculated and the 
atomic noise is generated from the probability distribu- 
tions of J y and J z . The error signal £(tk) is found and 
a frequency correction Auj(tk) is generated. The noise 
spectrum of the slaved oscillator, see Fig. lb, clearly 
shows that while for short time scales (< T, high fre- 
quencies) the noise is given by that of the free-running 
oscillator, at longer time scales (lower frequencies) the os- 
cillator is locked to the atoms and the remaining (white) 
noise is determined by the atomic fluctuations. The low- 
frequency white noise floor determines the long-term sta- 
bility of the clock and is the quantity we seek to optimize. 
In Fig. 1(c) we compare our analytical results with the 
results of the numerical simulations as a function of Ram- 
sey time T, and in Fig. 2(a) we show the scaling with 
the number of atoms. The analytical and numerical ap- 
proaches are in excellent agreement. 

So far we have assumed linear feedback and white 
noise; we now relax these assumptions. The stability 
limit identified above is mainly determined by the break- 
down of the assumption of small (i.e., linear) phase flue- 
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FIG. 2: Inverse fractional frequency stability l/o\, (arbitrary 
units) vs. number of atoms N, with Ramsey time optimized 
for a) white noise and b) 1// noise. Points: numerical simu- 
lations, lines: analytical results. Uncorrelated atoms (o) and 
optimal spin squeezed atoms (o) , both for linear feedback (full 
lines, open symbols) and nonlinear feedback (dashed lines, 
filled symbols). 



tuations. In fact, the stability can be improved consider- 
ably by using a feedback Auj which is a nonlinear func- 
tion of the error signal £. To investigate this we have 
included a nonlinear feedback Auj(tk) oc arcsin[£(tfc)/ J] 
in our numerical simulations. In Fig. 1(c) it is seen that 
nonlinear feedback performs better, and that it extends 
the validity of Eq. (JjjJ all the way to 7T ~ 0.1. For larger 
7T, the feedback loop fails, resulting in a rapid decrease 
in stability. If we optimize the Allan deviation in Eq. (JSJ 
for nonlinear feedback, under the condition 7T < 0.1, we 
find that the optimally squeezed states have A J y ~ iV 1 / 3 
(£ - TV -1 / 6 ) resulting in a stability scaling as N 
This represents again a relative improvement in scaling 
of TV 1 / 6 compared to the uncorrelated state for which the 
stability scales as iV" 1 / 2 . Detailed derivation of these re- 
sults will be presented elsewhere flj| . 

The assumption of white noise (<5</> 2 ) = 7T, is conve- 
nient for theoretical calculations, but in practice very- 
low-frequency noise is likely to have nontrivial spectrum 
such as 1// noise. To find the scaling with the number 
of atoms in this situation, we replace (6<fi 2 ) = 7T with 
the behaviour expected for 1// noise: (Sip 2 ) ~ (7T) 2 . 
Repeating all the calculations above we again find an im- 
provement by a factor of TV 1 / 6 by using squeezed states 
for the nonlinear feedback loop, and a factor of TV 5 / 24 
for linear feedback. In Fig. 2b we compare these scal- 
ing arguments to the numerical simulations and the two 
approaches are seen to be in very good agreement. 

To summarize, we have shown that entanglement can 
provide a significant gain in the frequency stability of an 
atomic clock when it is limited by the stability of the os- 
cillator used to interrogate the atoms. The optimal sta- 
bility is achieved by using moderately squeezed states, 
with a relative improvement that scales approximately 
as TV 1 / 6 with the number of atoms. These results are 
in contrast to previous studies using simplified decoher- 
ence models, which found that no practical improvement 
can be achieved with entangled states. Finally, we note 
a few interesting questions raised by our work. First, it 



would be interesting to see if there exists special quantum 
states of atoms and feedback mechanisms which optimize 
the performance of the clock. Second, the present results 
highlight that it is essential to have a realistic model of 
the noise (and possible stabilization mechanism) present 
in specific realizations of quantum information protocols. 
Although the protocol considered in this Letter exploits 
entanglement to stabilize a classical system (the local os- 
cillator), it would be interesting to study how similar 
considerations (e.g. 1// noise and collective decoher- 
ence) affect protocols such as quantum error correction 
codes , which use entanglement to stabilize a quantum 
system and protect it from decoherence. 
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